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Abstract 

We report on calculations of the heavy quark potential in SU (3) lattice gauge theory. 
Full SU(3) results are compared to three cases which involve gauge-fixing and projec- 
tion. All of these start from the maximal abelian gauge (MAG), in its simplest form. 
The first case is abelian projection to U(l) x U(l). The second keeps only the abelian 
fields of monopoles in the MAG. The third involves an additional gauge- fixing to the 
indirect maximal center gauge (IMCG), followed by center projection to Z(3). At one 
gauge fixing/configuration, the string tensions calculated from MAG U(l) x U(l), MAG 
monopoles, and IMCG Z(3) are all less than the full SU(3) string tension. The projected 
string tensions further decrease, by approximately 10%, when account is taken of gauge 
ambiguities. Comparison is made with corresponding results for SU{2). It is emphasized 
that the formulation of the MAG is more subtle for SU(3) than for SU(2), and that the 
low string tensions may be caused by the simple MAG form used. A generalized MAG 
for SU(3) is formulated. 
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1 Introduction 



One of the most important ideas in the theory of quark confinement is the notion of abelian 
projection, introduced long ago by 't Hooft 0. The basic conjecture made by 't Hooft was 
that in an appropriate partially fixed gauge, confinement can be understood for SU(N) 
pure gauge theory by retaining only the gauge fields belonging to the maximal abelian 
subgroup [/(l)^ -1 . In 't Hooft's original paper, there were several suggestions for how 
to partially fix the gauge, leaving the maximal abelian subgroup unfixed. Subsequent 
work in SU(2) lattice gauge theory has focussed on the well known 'maximal abelian 
gauge' (MAG). The main purpose of the present work is to perform lattice gauge theory 
calculations which test 't Hooft's idea of abelian projection for SU(3), using the simplest 
possible form of the MAG. 

1.1 Beyond Abelian Projection 

The projection process can be carried beyond abelian projection, keeping a still smaller 
set of degrees of freedom. Closely related to abelian projection or abelian dominance is 
the postulate of monopole dominance. This assumes that the confining physics in the 
abelian gauge fields originates in the magnetic current of monopoles. The monopoles are 
assumed to be similar in nature to 't Hooft-Polyakov monopoles, moving on world lines. 
It is known that a dense network of magnetic current is the agent of confinement in U(l) 
lattice gauge theory f2|. The key question is whether a corresponding network can explain 
confinement in SU (N) lattice gauge theory. 

Another possibility is that the confining degrees of freedom are really associated with 
the Z(N) center of SU(N). Since Z(N) C U(1) N ~ 1 , this can be studied after making an 
abelian projection to U(1) N ~ 1 . After the initial gauge-fixing and abelian projection, the 
aDe ij ail links are subjected to a second gauge-fixing and projection, leaving a set 
of Z(N) or center- valued links. This is the 'indirect maximal center gauge' (IMCG). The 
IMCG uses the MAG as an intermediate stage, as opposed to the 'direct maximal center 
gauge' (DMCG), which does not. For either case, center dominance is tested by using 
Z(N) links, or closed world sheets of plaquettes carrying Z(N) flux. The latter are called 
(P)-vortices. 

1.2 Gauge Ambiguities and SU(2) vs. SU(3) 

In this paper, we work exclusively with 'maximal' gauges; gauges in which a certain 
functional of the links is maximized before the links are abelian or center-projected. The 
use of a maximal gauge has an intuitive appeal. In the case of abelian projection, if the 
charged degrees of freedom are to be deleted, why not first transform to a gauge where 
they are made as small as possible? Maximal gauges also have the desirable practical 
feature that the noise level in projected configurations is quite low. 

A complication associated with maximal gauges is the presence of gauge or Gribov 
ambiguities. The functionals used have many local maxima for a given configuration. Each 
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local maximum of the gauge functional is a 'gauge copy', and the results of numerical 
calculations depend on which gauge copy is chosen for each configuration. That this 
variation with gauge copies occurs at all is caused by projection. If the complete set of 
gauge degrees of freedom is kept, a gauge invariant quantity like a Wilson loop has the 
same value at any point on the gauge orbit of a configuration. But if certain degrees 
of freedom are deleted, then the value of the projected Wilson loop will in general be 
different at different points on the orbit, even if the chosen points are local maxima of a 
gauge functional. 

SU(2) calculations have shown a clear qualitative trend for this variation due to gauge 
copies. Gauge copies with higher values of the functional produce projected configurations 
which are more ordered and produce lower string tensions. This is true for the maximal 
abelian gauge followed by abelian projection, as well as the direct maximal center gauge 
followed by center projection. In both cases, if for each configuration a randomly chosen 
gauge copy is used, the projected string tension turns out to be larger than the full 577(2) 
result. If instead, a collection of gauge copies is generated for each configuration, selecting 
a subset and then picking the member of the subset with the highest value of the gauge 
functional leads to a reduced string tension. For the MAG and abelian projection, recent 
work suggests that use of the 'best' (global maximum) gauge copy for each configuration 
would lead to an abelian string tension quite close to the full SU(2) value 0. For the 
DMCG and center projection, there is more variation. As gauge copies with increasing 
values of the center gauge functional are used, the Z(2) string tension decreases from 
an initial value above the full SU(2) string tension to a value well below it |4]]. While 
the global maximum of the gauge functional evidently plays a different role in the two 
projection schemes, for the purposes of this paper, the important point is that for both 
abelian and center projection, there exist gauge copies that give string tensions very close 
to the exact SU(2) value. 

For SU(3), the correct formulation of the MAG condition is a more subtle matter than 
it is for SU(2) . This is discussed in detail in Sec.(0). The fact that there is more than 
one possible form for the MAG is important, since one of our main results is that for the 
simplest form of the MAG in SU(3) , it is practically certain that no gauge copies exist 
which produce a U{1) x U{1) string tension close to the full SU{3) value. 

Maximal center gauge conditions are discussed in Sec.(||]). In our numerical work, we 
used the IMCG which has the MAG as a preliminary step. Here we also find a low string 
tension. The amount by which the Z{3) string tension is low depends on the treatment 
of gauge copies. Since the MAG and abelian projection are an intermediate stage, the 
low values we find for the Z{3) string tension are not necessarily a problem of center 
projection in SU(3), but may be reflecting the difficulty discussed above for the simple 
form of the MAG. 

Our numerical calculations were performed at /3 — 5.9 on a 10 3 x 16 lattice and 
at /3 = 6.0 on a 16 4 lattice. The full SU(3) results are in Sec.(|4]). For the MAG, the 
U{\) x U{\) results are in Sec.(|), and monopole results are in Sec.(|]). For the IMCG, 
the Z{3) results are in Sec.([7|). We find that as in SU(2), the monopole and center vortex 
degrees of freedom are intimately connected. At f3 — 6.0, for both the MAG and IMCG, 
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we investigated the dependence on gauge copies. We confirm for both gauges, that gauge 
copies with higher values of the gauge functional lead to lower string tensions. This is the 
same general behavior as seen in £77(2). Sec.(|8|) contains a summary and conclusions. 

2 Maximal Abelian Gauge 

In this section, we discuss the maximal abelian gauge condition for SU(3), comparing to 
the familiar case of SU(2). For £77(2), we regard as the abelian field, and by making 
gauge transformations, suppress the 'charged' fields A 1 and A 2 ^ to the greatest extent 
possible. The continuum functional 

G™1 = J d"x[{Alf + [Alf] (1) 

is minimized over £77(2) gauge transformations. Near a minimum, G™% g is stationary 
when the gauge fields are varied by an infinitesimal gauge transformation, 

Al -> Al - d^ a + e abc 9 b A c ^. (2) 

Since G^ 2 g is invariant to a small abelian gauge transformation involving # 3 , only trans- 
formations involving O 1 ^ 2 lead to gauge conditions. If we combine the gauge fields into 
charged combinations, A^ = (A* ±iA^)/\/2, the conditions can be written as 

d„A+ ± iAlA± = 0. (3) 

Turning to ££7(3), there are now two abelian gauge fields, A^ and A* . The remaining 
six gauge fields are all charged with respect to at least one of A^ , A % . Treating these 
charged fields democratically, the simplest maximal abelian gauge condition for ££7(3) 
would minimize the functional 

Gmag = J d'xUlf + {Alf + (4) 2 + (4) 2 + (Al) 2 + (Al) 2 } (4) 
over the £77(3) gauge group. An infinitesimal ££7(3) gauge transformation is 

ai^ai- o,e a + r bc e h A^ (5) 

where f abc are the ££7(3) structure constants. Since G^ g is invariant to a small abelian 
gauge transformation involving 8 3 and 8 8 , there are six gauge conditions arising from the 
remaining six angles. These conditions can easily be expressed in terms of the three ££7(2) 
subgroups of ££7 (3) . We denote the three subgroups and their gauge fields as follows: 

I- Al,Al,Al 

V '■ ^% 
11 ■ 4 6 A 7 A 3 " 
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where A^' = + ^4®, Af" = -\A^ + ^A^ } and with some abuse of notation, quark 



su3 
mag 
su3 
mag 



model terminology has been used to denote the subgroups. Consider the response of G 
to an infinitesimal gauge transformation in the /-spin subgroup. Only the terms in G 
involving A 1 , A^ vary under this transformation, so demanding that G^ g be stationary 
leads to Eq.© again. Likewise, demanding that G^ g be stationary with respect to a 
V^-spin gauge transformation leads to 

V£'±iAj'Aj' = 0, (6) 

where At' = (A* ± iA^)/y2. Finally, [/-spin gauge transformations lead to 

9 M 4"±i4"4" = 0, (7) 

where A^" = (A 6 ^ ± iAjJ/y/2. So the SU(3) maximal abelian gauge functional Eq.([|) 
leads to the requirement that the gauge fields be 'maximally abelian' with respect to each 
of the three SU{2) subgroups of SU{3). 

Minimizing the continuum functional Eq.([|) is equivalent to maximizing the lattice 
functional G s ^ g defined by 



2E(l(^)ll| 2 +K^)22| 2 +l(^ 



|2 



M/331 



X.fl 



We will use Eq.([8D for our MAG gauge- fixing in this paper. This same functional has 
been used in the small number of other papers published on the MAG for SU(3) [0. [|. 

At a local maximum of Eq.(^), three conditions are satisfied. These are the lattice 
analogs of the continuum conditions Eqs.(Q),(|6]), and (|7]), and are satisfied when G^ g 
is stationary with respect to gauge transformations in the /, U, and l^-spin subgroups. 
To write the conditions explicitly, we define three site-valued matrices. For the /-spin 
subgroup, the matrix X(x) is 

X(x) = E [u»(x)\ 3 Ul(x) + U\{x - fiXJJ^x - £)] • (9) 

For the V and t/-spin subgroups, the corresponding matrices X' and X" are obtained by 
replacing A 3 in Eq.@ by A' 3 and A3, respectively, where 

K = \h + ^A 8 , AS = -~A 3 + ^A 8 . (10) 

When G™a g is stationary against gauge transformations in the /-spin subgroup, the 12 
and 21 matrix elements of X(rr) must vanish; similarly, for the \^-spin subgroup, the 
13 and 31 matrix elements of X'(x) must vanish; and for the £7-spin subgroup, the 23 
and 32 matrix elements of X"(a;) must vanish. In practice, we require that these matrix 
elements be very small. The numerical methods we use to reach these conditions with 
high accuracy will be discussed in Sec(^). 
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2.1 Generalized Maximal Abelian Gauge 



So far there appears to be a close analogy between the MAG conditions for SU(2) and 
SU(3). However, there are more general possibilities for SU(3). To see them, we first 
note that G™% g is in the form of a mass term for the gauge fields Ai, A 2 . Such a term 
would arise in a gauge-Higgs theory, with an adjoint Higgs field $ a coupled to the SU (2) 
gauge field. If the Higgs field has a vacuum expected value along the 3-axis, equal mass 
terms for A 1 ^ and A 2 ^ will be generated. Take the nondynamical Higgs field to have a fixed 
length, $ a $ a = $q, and write the continuum Higgs action 



n 



su2 
higgs 



I I d 4 xD„$ ■ £>„$, 



where (D^) a = d^ a + e abc A b ^ c . In the Higgs gauge with $ = $ 3, becomes 

(K) 2 +K) 2 ) d ' x ( 12 ) 



risul 
higgs 



*8 



which is up to a constant, the same as G^ 2 g . 

The SU(3) functional Eq.(f|) is also in the form of mass terms for the charged gauge 
fields. It is of interest to see how the functional Eq.([|) is related to the action of a fixed 
length adjoint Higgs field $, with components $ a . We write the continuum action 



/^isu3 
higgs 



d xD„<$> ■ DJb, 



(13) 



where (_D M $) a = <9 M $ a + f abc A b ^ c . For SU(3), going to a Higgs gauge means putting 
the Higgs field in the Cartan algebra, $ = $o(3cosx + 8sinx), or in matrix form 
$ = $0(^3 cos X + -^8 sm x)/2- The value of the angle x cannot be changed by an SU (3) gauge 
transformation, so there are distinct 'moduli' or Higgs gauges parameterized by \. This 
is a difference between SU(2) and SU(3) . In a generalized Higgs gauge where we allow x 
to depend on x, the functional Eq.flT3"|) takes the form 



f~tsu3 
higgs 



$ 2 / d*x{d, X d»X+ cos^(x) (At) + (A 



cos 2 (^ - x) [(KY + (4) 2 ] + cos 2 (^ + x) [(Alf +(Al) 2 ] } • (14) 



Even for x independent of x, this is clearly different from Eq.(^). No individual adjoint 
Higgs field can give equal coefficients for the various charged vector fields in the func- 
tional. Formally, equal coefficients for the gauge fields can be obtained by taking x to be 
independent of x, and then averaging over x i n Eq. flljj) . 

Minimizing the functional Eq.flT4|) in the continuum is equivalent on the lattice to 
maximizing the functional Gf^g gs defined by 



G higgs = H tr (( A 3Cosx + A 8 sinx)t/ M (a:) (A 3 cosx + A 8 sin x'Wfa)) 



(15) 



X,fl 
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where \ — x{ x ) an d x' = x{ x + A a )- If X is taken to be independent of x, and Gf™g gs 
is averaged over x, ^higgs reduces to G^ g . Although not used in the present paper, the 
functional G^L would appear to be an interesting alternative to G^f g , particularly for 
the case where \ is allowed to depend on x. We intend to explore this in future work. 



2.2 Abelian Projection 

After the MAG gauge-fixing process is complete on a given configuration, each link 
Up is factored into a diagonal £7(1) x £7(1) part u^, times a 'charged' part U^ h , so 
Ufj, = uJJ^ 1 . The has generators in the Cartan algebra and represents the 'phase' 
of the full SU(3) link matrix U^. Short-range physics is hopefully isolated in the matrix 
U^ h , which depends on the remaining six angles. In contrast to the SU(2) case, simply 
extracting the phases of the diagonal elements of £7^ will not result in a matrix with 
det(w M ) = 1. In this paper, we use two different methods to extract w M . In the 'symmetric' 
method, the phase of a given element of is set equal to the phase of the corresponding 
diagonal element of U^, minus the average of the three diagonal phases. In the second, 
'optimal' method, we find the £7(1) x £7(1) link which maximizes Re(tr(u M t/r) |7j. This 
finds the 'best' £7(1) x £7(1) approximation to £7^, and requires an iterative procedure. 
We have carefully checked that results produced by symmetric or optimal methods differ 
by a negligible amount, considerably smaller than statistical errors. 

Once the £7(1) x £7(1) links w M are determined, the abelian projection is carried out 
by setting U^ h to the identity on every link of the lattice, so the resulting SU(3) links 
are equal to the u M . Using the in place of the directly tests abelian dominance for 
S77(3). The results of these calculations are reported in Sec. (|5|) . 



2.3 SU(3) Monopoles 

In this section, we discuss the extraction of magnetic currents and calculations with 
monopoles for SU(3), after reviewing the SU(2) case ||. There, after MAG gauge-fixing 
and abelian projection, the links are diagonal 2x2 matrices w M given by 



u, 








(16) 



where (j>jf 



1,2 and 



+ 



e (-7T,+7T), J _ ^ 

procedure B to plaquettes formed from the 



However, m l + 



m. 



^ = 0- Applying the Toussaint-DeGrand 
bjj gives two integer-valued magnetic currents 
is satisfied explicitly on every link, so there is only one 



magnetic current for SU(2). 

For SU(3) , after going to the MAG and projecting to £7(1) x £7(1), the abelian links 
are of the form 

/ exp(z^) \ 

exp(z^ 2 ) 
V exp(z^ 3 ) J 



u, 



(17) 
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with 4>]} + 4>ff + = 0, mod(27r). Our procedure for extracting magnetic currents closely 
follows that just discussed for SU(2) . Plaquette angles are computed for each color, and 
the Toussaint-De Grand procedure is applied to each. This produces a magnetic current 
for each color, m^, j = 1,2,3. Because of the mod(27r) in the sum of the angles <fiff , 
the total magnetic current + + does not vanish link by link. However, as 
shown in Sec.(^), the total magnetic current has no physical effect, so there are actually 
only two independent physical magnetic currents. This method of treating the magnetic 
current has the advantage that it does not force a particular subgroup structure for SU (3) 
monopoles. Evidence that they are in fact associated with £77(2) subgroups is presented 
elsewhere [ |IO|j . 

The monopole Wilson loops are calculated exactly as in SU(2) or U{1) @. Each color 
of magnetic current produces a lattice magnetic vector potential A 3 m , and field strength 
F^, = d^A 3 mv — d+A J m . The monopole Wilson loop is then obtained from the flux of 
the dual of through the loop. We have 

{WLon) = (exp (*27T £ D^x)*F^(x)\ ) , (18) 

where is unity on any simple surface of plaquettes whose edge is the Wilson loop, and 
< > m means the average over the ensemble of magnetic currents. The dual field strength 
is defined by F^^ = e ilvol ^Fl nOL ^j2. Before extraction of potentials, the < > are 

averaged over color. 



3 Maximal Center Gauge 

Gauge transformations involving a discrete subgroup are well-defined when acting on a 
path ordered exponential of gauge fields, 

U(xt, x 2 ) = P (exp(J A^dx^J . 

Hence a gauge can be sought where a set of U{x\, x 2 ) are as close as possible to elements 
of a discrete subgroup of the gauge group. 

In lattice gauge theory, the U{x\,xz) are taken to be the links of the lattice. The direct 
maximal center gauge for SU(N) lattice gauge theory seeks maxima of the functional 

Gs£, = El tr ^l a . ( 19 ) 

x,n 

over the SU(N) gauge group. Since trL^, is real for 577(2), an alternate form for SU(2) is 

GTl 9 = E i^U,) 2 . (20) 
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The Up are in the fundamental representation in all the formulae above. Fundamental 
and adjoint traces are related by \r(U^) = |trL^| 2 — 1, so Eq.([L9|) can be rewritten as 

Gi = EK 4 ) + l). (21) 

X,fJ, 

which shows that the DMCG can be thought of as an adjoint Landau gauge. Gauge 
transformations are made so as to make the gauge-fields as small as possible, but any 
Z(N) factors in the fundamental links are left unfixed. 

Although we will restrict ourselves to the bilinear condition Eq.([L9[) in this work, it 
is worth noting that for N > 2, there are additional Z(N) invariant functionals of the 
links that can be constructed |TT[] . For example in SU(3) , the expressions tr (U^U^U^), 
tr(6 r M 6 r /i ) • trf/^, and ^(6^) • tr([/ M ) • tr({7 M ) are all Z{3) invariant. From these a family 
of real-valued, Z(3) invariant functionals of the links could be written down and used to 
define 577(3) maximal center gauges. It would be interesting to know if such gauges differ 
in their predictions from Eq. flTP|) . 



3.1 Indirect Maximal Center Gauge 

In this paper, we take the indirect route to the maximal center gauge for SU(3). This 
is the same method first used for Z(2) center projection in SU(2) |12j|, and fits naturally 
with our study of the MAG for SU(3). A desirable feature of the IMCG is that it allows a 
study of the relationship between monopoles and center vortices. We start by going to the 
MAG as discussed in Sec.(@), using the functional of Eq.(||). The resulting 577(3) links 
are then projected to U{\) x £7(1) links using the methods discussed in Sec.(|2.2|). 
Parameterizing the w M as in Eq. (|T7D , and using them to evaluate Eq. (|T9D , the functional 
to be maximized for the IMCG becomes 

£ [3 + 2cos(0i 1 (x)-0 22 (x)) + 2cos(0i 1 (x) + 20^ (x)) (22 ) 

X,fl 

+2cos(0 22 (x) + 20i 1 (x)); , 

where the unitarity condition has been used to eliminate 33 . The gauge functional of 
Eq.(|2"2"|) is maximized by successively applying abelian gauge transformations at each site. 
The gauge transformations are parameterized as 

where 

A 8 = diag(l, 1, -2), A 3 = diag(l, -1, 0). 

Maximizing at a site x gives two coupled equations in a 3 (x) and a$(x). These equations 
are not algebraically soluble, and their solution is found iteratively. Making the required 
gauge transformation at x generally causes the maximization conditions at the nearest 
neighbors of x to be violated, so reaching a local maximum of Gf^ requires repeated 



f~i su3 



8 



cycles through the lattice. At a maximum, two gauge conditions are satisfied. The first 
comes from the Ag term in the U(l) x U(l) gauge transformation, 



£(9; [sin(20i 1 (x) + <pf(x)) + sin(20j 2 (x) + <ff(x)j\ = 0. (23) 



i 1 

The second comes from the Ac. term, 



E^sin^x) - <j>f{x)) + sin(20i 1 (x) + <t>f{x)) - sin(2# 2 (s) + <ft(x))\ = 0. (24) 



In the naive continuum limit these two conditions reduce to the abelian Landau gauge 
conditions 

d,Al = 0, d,Al = 0, 

where aA 8 ^ = {<f)j}{x) + <f) 22 (x))/2, and a A 3 = {(j^Hx) — <j) 22 (x))/2. In our lattice calcula- 
tions, the right hand sides of Eqs.(|23|) and ( f2"4"|) are not actually zero, but are very small. 
This is discussed in Sec.(^). 

After a local maximum is reached to sufficient accuracy, the links are projected to 
the nearest element of the center group, Z(3), where 



Z{3) G {exp(-i27r/3),l,exp(+i27r/3)} 



x 1. 



This results in a set of Z(3)-valued links, z^(x). These are used to compute the Z(3) 
Wilson loops discussed in Sec.([7|). 

Since our results for center projection are all for the IMCG, it is worth commenting on 
possible differences between the direct and indirect maximal center gauges. The IMCG 
depends on the specific functional form used for the MAG functional, which as discussed 
in Sec.(0), is not unique for SU(3). This means the IMCG results obtained in this paper 
using the Eq.(^|) form of the MAG are not necessarily indicative of those that would be 
obtained in the DMCG for SU (3) . For the case of an SU (2) gauge group, where the MAG 
functional is unambiguous, we have found no statistically significant difference between 



indirect and direct methods of reaching the maximal center gauge [pfl . A similar result 
may hold in the end for SU(3), but this will probably require an MAG functional different 
from Eq.(||). 



4 SU(3) Results 

Our calculations were performed at couplings (3 = 5.9 and 6.0 on lattices of size 10 3 x 16 
and 16 4 , respectively. The correlation length as determined by the string tension is ap- 
proximately 4a at (3 = 5.90 and 4.5a at (3 = 6.0, so at (3 = 5.90, our 10 3 x 16 lattice is over 
twice the correlation length in size, while at (3 = 6.0, the 16 4 lattice is more than three 
times as large as the correlation length. Finite size effects should be small, particularly at 
(3 = 6.0. In each case, SU(3) configurations were generated using the Wilson action with 



a heatbath update and overrelaxation |I4"||. The first several thousand configurations were 
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discarded to allow for equilibration, and every 20 th configuration after that was used for 
measurements. We analyzed 440 configurations at f3 — 5.9, and 600 at (3 = 6.0. 

In this section, we present full SU (3) results for the heavy quark potential. Results for 
the MAG U(l) x U(l) , MAG monopole and IMCG Z(3) potentials are presented in the 
sections immediately following. The full SU{3) potential is of course gauge- invariant, and 
therefore not subject to the uncertainties of gauge-fixing and projection. On the other 
hand, the noise level of SU(3) Wilson loops is much greater than for projected loops, and 
variance reduction is essential to extract the potential. The presence of a large Coulomb 
term also complicates the determination of the string tension. For the SU (3) potential, 
modified Wilson loops were computed using a multihit variance reduction scheme [13] and 
spatial smearing ]TJ|. In computing the heavy quark potential, we used loops of spatial 
extent R = 2a to R = 5a {(3 = 5.9), R = 2a to R = 7a (/3 = 6.0). The maximum temporal 
extent was T = 12a in both cases. From the Wilson loops, we extracted potentials by 
fitting — \n(W(R, T)) to a straight line in T, the slope giving V(R) (T^j. The values of 
T included in the fits for each R were T = R + a to T — 12a. The potentials were then 
fitted to the form V(R) = V + a/R + aR. 

Our results for the full SU(3) potential parameters are presented in Table [D, For 
the purposes of this paper, the string tensions in Table [I] are the most significant, since 
the various projection schemes we will compare to in later sections do not (nor are they 
expected to) reproduce the Coulomb part of the potential. Comparing our SU(3) results 
to the modern high statistics calculations available in the literature, our (3 = 6.0 string 



tension is within error bars of the values quoted in [Tj|, [L9], |2(J, j21| . At (3 = 5.90, the 
only calculation with comparable accuracy is p2] , where the value reported is somewhat 
higher than ours, by an amount slightly outside of error bars. However, it should be noted 
that at (3 = 6.0, the string tension quoted in |22| is also somewhat higher than the results 
in [P| [T9|, EjOL EHJ. We conclude that our full SU(3) string tensions are of good quality, 



comparable to the best values available in the literature. 



5 U{1) x U{1) Results 



This section contains our most important results. Abelian projection after gauge-fixing to 
the MAG retains the largest number of degrees of freedom of the various projections we 
consider, and so is a less severe truncation of the theory than retaining only the abelian 
fields of monopoles, or keeping only Z(3)-valued links after center projection. 

The process of actually finding a local maximum of the MAG functional was car- 
ried out iteratively, determining the gauge transformation g(x) at a specific site which 
maximizes the terms in G^ g that are affected by g(x). Gauge transformations were 
successively chosen from the three independent subgroups of SU(3). An arbitrary SU(3) 
gauge transformation can be expressed as a product of such SU (2) gauge transformations. 
To improve convergence, we used an overrelaxation procedure, replacing g(x) with g u} (x) 



2vJ. For gauge transformations which are elements of SU(2) subgroups, it is possible to 
work out g u3 {x) exactly, rather than use the Taylor expansion approximation described in 
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It was found that uo = 1.85 was the optimum overrelaxation parameter when using 
the exact expression for g UJ (x). Referring back to the conditions discussed after Eq.(|8]), 
gauge fixing iterations were performed until 

(|X(x) 12 | 2 + |X(x) 21 | 2 + |X'(x) 13 | 2 + |X'(x)3i| 2 + |X''(x) 23 | 2 + |X''(x)3 2 |^ < KT 12 , (25) 

where <> implies a lattice average. Eq.fl2"5|) is a very tight constraint. We are confident 
that when it is satisfied, the configuration is very close to a local maximum of G^f g . As 
with any maximal gauge, the functional G^ g has many local maxima or gauge copies. 
A single gauge-fixing brings a configuration near a random local maximum of the func- 
tional, unrelated in general to the global maximum. Actually finding the global maximum 
is impractical. However, it is feasible to generate a sample of gauge copies for each con- 
figuration, and study how the projected string tension varies when copies with steadily 
increasing values of the functional are used. Our results are presented first for the case of 
one gauge copy/configuration, followed by analysis of the effect of gauge ambiguities. 

After reaching a local maximum for each configuration, the SU(3) links were projected 
to the U(l) x £7(1) links u^, using the methods discussed in Sec.( [2.2|) . Spatial links were 



smeared using the Ape method [16] with a staple weighting parameter of 0.7, iterating 



the procedure 8 times per link. Then abelian Wilson loops were computed using 

(26) 



These loops were used to compute the U(l) x U{1) heavy quark potential V(R) at both 
(3 = 5.90 and (3 = 6.0. The results for fits to V(R) of the same type as for the full SU(3) 
potential are shown in Table [| 

The U(l) x U(l) string tensions shown in Table H are for one local maximum or gauge 
copy /configuration, and are smaller than the corresponding full SU(3) results, by 10% 
at (3 = 6.0. This is in contrast to the behavior of SU(2) in the MAG for one gauge 
copy /configuration. There similar gauge-fixing methods produce U(l) string tensions at 
least 10% larger than full SU{2) at (3 SU{2) = 2.4,2.5,2.6 |g g|. 

Turning to gauge ambiguities, in order for their effect to bring the U(l) x U(l) string 
tension closer to the full SU(3) answer, local maxima of the gauge functional with higher 
values of the functional would have to produce higher abelian projected string tensions. 
This is the opposite of what happens in SU(2). To explore whether this unexpected 
behavior occurs for the SU(3) MAG functional of Eq.(f|), we generated N copy = 10 gauge 
copies/configuration for each of 200 configurations at (3 = 6.0. To generate a gauge copy, 
a move is made to an arbitrary point on the gauge orbit of a configuration by making a 
random gauge transformation. The gauge-fixing algorithm then takes the configuration 
to a nearby local maximum. Due to abelian projection, the U(l) x U(l) string tension 
is dependent on which copies are used for each configuration. Our treatment of the 
variation due to gauge copies follows the method of Bali, et. al. |24|], developed for the 



SU(2) case. For each configuration, a subset of size n g of the N copy copies is chosen. Only 
the configuration in the subset with the highest value of the gauge functional is retained 
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and used to calculate Wilson loops. The Wilson loops are averaged over all such subsets 
of size n g . For n g — 1, this is just the average over all N copy copies for each configuration, 
whereas for n g = N copy , only the copy with the highest value of the functional for all N copy 
copies is used for each configuration. The resulting Wilson loops are then averaged over 
configurations and used to calculate potentials. The output of this analysis is a set of 
values of the £7(1) x £7(1) string tension for n g ranging from 1 to N copy . The results are 
shown in Table |[ As can be seen there, the U(l) x £7(1) string tension steadily decreases 
as n g increases. At n g = N copy = 10, the U{1) x £7(1) string tension has decreased by 
10%, and the ratio of the £7(1) x 17(1) to the £77(3) string tension is 0.80(4). 

We have thus found that the effect of gauge ambiguities in £77(3) is similar to SU(2), 
namely, the higher the value of the MAG functional, the lower the abelian projected string 
tension. In SU(2), this trend is beneficial, since for a randomly chosen gauge copy, the 
£7 (1) string tension is above the full SU (2) result, and taking gauge copies with successively 
higher values of the functional brings the string tension down, closer and closer to the 
full ££7(2) value. For SU(3), using the Eq.(§ form of the MAG, the £7(1) x £7(1) string 
tension from a random gauge copy is too low to start with, and gauge copies with higher 
functional values produce U(l) x £7(1) string tensions still further below the SU(3) string 
tension. 



6 Monopole Results 

Before presenting results on the heavy quark potential due to monopoles, we discuss a 
point mentioned in Sec.( [2.3| ), related to the determination of the magnetic current. Our 
procedure for locating monopoles produces a separate current for each color, j = 1, 2, 3. 
This method treats the three colors symmetrically, and is not prejudiced toward a specific 
group structure for the monopoles. Only two of these currents are really independent. 
Although it does not vanish link by link, the total current, + + m^, is in effect 
a null current. We have checked this by using the total current to calculate monopole 
Wilson loops. The resulting string tension (at (3 = 5.90) is zero to four decimal places, 
a^ on = 0.0000(4), as it should be. 

Turning to the monopole potentials, at (3 — 5.9 and 6.0, we computed monopole Wil- 



son loops for each color using Eq. (p~8|) , averaged the results over color, and extracted 
heavy quark potentials using the same fit methods discussed earlier. The results are 
presented in Table [| along with the current density, or fraction f m of links carrying mag- 
netic current. The results shown in Table |] are for one MAG gauge-fixing/configuration, 
and are on the same footing as the £7(1) x £7(1) results in Table 0. The inequality 
&su(3) > &u(i)xu(i) > &mono is seen to hold for both (3 values, the monopole string ten- 
sions being approximately 25% smaller than the SU(3) string tensions, a discrepancy well 
outside of error bars. 

It is of interest to compare these results with corresponding ones for £77(2). The full 
SU(2) string tensions at Psu(2) = 2.4 and (3su(2) = 2.45 are numerically close to the full 
£77(3) string tensions at (3su(3) = 5.9 and f3su(3) = 6.0, respectively, implying similar 
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correlation lengths. However, for 577(2), the monopole string tensions for one gauge- 
fixing/configuration are within 5% of the full SU(2) values ||. Further, the magnetic 
current densities f m , are more than twice as large as the corresponding SU(3) results 
shown in Table (| . Thus for SU(3), the MAG condition of Eq.(|8]) leads to projected 
configurations with a severe shortage of monopoles. Despite this shortage, the magnetic 
current density shows scaling properties rather similar to those found in SU(2). If the 
magnetic current density is a physical density, dividing it by the cube of a physical energy 
should produce a scaling quantity, i.e. one independent of (3. Computing / m /(o" mono ) 3 / 2 
for (3 = 5.90 and (3 = 6.0 gives 1.1(1) and 1.0(1), respectively. Corresponding results 
for SU(2) at (3 SU {2 = 2.40 and (3 SU{2 = 2.45 are 1.6(1) and 1.7(1), respectively §. So 
although the magnetic current density for £77(3) is significantly smaller, it scales with f3 
in a similar way to the SU(2) case when the monopole string tension defines the unit of 
energy. 

Just as they affect the U{\) x £7(1) potential, MAG gauge ambiguities affect the 
monopole potential parameters. Using the same set of N copy = 10 gauge copies/configuration 
at (3 — 6.0, we can determine how the monopole string tension and magnetic current 
density vary with the number of gauge copies. The analysis proceeds exactly as in the 
U{\) x £7(1) case. The results are presented in Table §. As n g increases, the monopole 
string tension decreases by essentially the same factor as the U(l) x £7(1) string tension, 
and the ratio of the magnetic current density to (a mon ) 3 ^ 2 remains approximately con- 
stant. By the time n g = N copy = 10 has been reached, the f3 = 6.0 value of the monopole 
string tension is 34% smaller than the S77(3) string tension. 

7 Z(3) Results 

To transform an SU(3) configuration to the indirect maximal center gauge requires two 
gauge-fixings and one projection. First a local maximum of the MAG functional is found, 
followed by abelian projection. The resulting U{1) x U (1) links are then gauge-fixed by 
finding a local maximum of the IMCG functional Eq.(^), using the procedure described 
in Sec. (3.1). Iterations were performed until the sum of the squares of Eqs.(p3|) and ( p4"D 
was less than 10~ 14 , when averaged over the lattice. During the gauge fixing we monitored 
the modulus of the abelian trace, |tra^|. If all links were elements of the center, the value 
of |(|tra M |) would be unity. At (3 = 5.90, we observed it to change from an average value 
of 0.525 before gauge fixing to a value of 0.922 after gauge fixing. The behavior at (3 = 6.0 
is similar. 

The IMCG has gauge ambiguities of two types; those inherent in the MAG, as well 
as those that are associated with the second gauge- fixing just described. We first present 
results for one gauge-fixing of each type/ configuration. This puts each configuration 
in a random local maximum of both gauge conditions. The resulting center-projected 
links were used to compute Z(3) Wilson loops and extract a heavy quark potential as 
in the previous cases. Similar to the monopole case, the noise level is very low, and the 
potential is almost perfectly linear. Table ||] shows the potential parameters, along with p, 
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the fraction of plaquettes with non-trivial Z{3) flux. These plaquettes are conjectured to 
represent the geometrical centers of finite-sized center vortices. Usually called P-vortices, 
we simply refer to non-trivial Z(3) plaquettes as vortices. From Table ^| and Table [I], the 
ratio of the Z{3) string tension to the full SU{3) value is 0.88(8) at f3 = 5.9, but falls to 
0.80(8) at (3 = 6.0. 

The vortex density is an areal density, so if it is physical, dividing by the square of 
a physical energy should produce a scaling quantity. The natural quantity to divide by 
is the Z(3) string tension. Computing the ratio of the vortex density to the Z(3) string 
tension gives .73(1) and .75(1) for (3 = 5.90 and (3 = 6.0 respectively. For the DMCG 
in SU(2) the results at Psu{2 — 2.40 and /3su(2 = 2.50 are quite similar; 0.76(1) and 
0.78(1), respectively Thus for both SU(3) and SU(2), dividing the vortex density by 



the Z(N) string tension removes most of the /9-dependence and produces an approximate 
scaling quantity. 

We also explored the effect of gauge ambiguities on the IMCG results. To avoid 
an expensive simultaneous treatment of both MAG and IMCG gauge ambiguities, we 
generated 10 IMCG gauge copies only for the MAG copy with the highest value of the 
MAG functional. This is the best MAG copy from the viewpoint of the MAG functional. 
The analysis of the Z{3) string tension proceeded in exactly the same manner as the 
U(l) x U(l) case. The results are shown in Table |7|, where it is seen that the Z{3) string 
tension decreases by O(10%) in going from n g = 2 to n g = 10. The ratio of vortex density 
to Z{3) string tension is approximately constant. At n g = 10, value of the Z(3) string 
tension is roughly 40% smaller than the (3 = 6.0 SU(3) string tension. 

If a different choice of MAG copy had been made, it is reasonable to expect a similar 
O(10%) variation due to IMCG gauge ambiguities. If for example, a random MAG gauge 
copy is used for each configuration and IMCG copies generated, the Z{3) string tension 
would be expected to decrease from its value given in Table |^ to one O(10%) lower, so for 
n g = 10, the Z{3) string tension would be O(70%) of the (3 = 6.0 SU(3) string tension. 

A reasonable conclusion is that using the gauge functionals Eq.(||) for the MAG and 
Eq.(|22|) for the IMCG, whatever choice is made of MAG or IMCG gauge copies, the result 
is a (3 — 6.0 Z(3) string tension whose central value is at least 20% smaller than the full 
SU (3) string tension. 

Finally, it is of interest to study the connection between vortices and monopoles. 
Suppose the magnetic current is non-zero on a dual lattice link. This link is dual to a 
3-cube on the original lattice. The partial magnetic current density f m {n) is defined as 
the current density for links where the dual cube has n faces pierced by Z(3) vortices. 
The total magnetic current density is then 



fm fn 



n 



n=0 



Since Z(3) center vortices form closed surfaces, / m (l) = 0, but all other possibilities can 
and do occur. The study was carried out for one gauge-fixing of each type/configuration. 
We computed the ratio r m {n) = f m (n)/ f m , which measures the fractional occurrence for 
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the cube dual to a magnetic current link to have n faces pierced by Z(3) vortices. The 
results are shown in Table §, and hold independent of color. Thus approximately 85% 
(83%) of the time at (3 = 5.9 (6.0), the cube dual to a magnetic current link has at least 
two faces pierced by vortices. This says there is an intimate connection between center 
vortices and monopoles in SU(3) , similar to what was found previously in SU{2) [13|, [25 |. 
It could in fact be true that the actual fractions are 100%, with the smaller numerical 
results caused by inaccuracies in locating monopoles and vortices on the lattice. 



8 Summary and Conclusions 

We have calculated the heavy quark potential for SU{3) lattice gauge theory at (3 = 5.90 
and f3 = 6.0. These SU(3) results have been compared with three cases which project onto 
a reduced set of degrees of freedom after gauge-fixing. The U{1) x U{1) and monopole 
results involve the version of the MAG expressed in Eq.(|8]), while the Z(3) results involve 
in addition the IMCG gauge- fixing of Eq.(p2|). In Figs.(|]) and (||), we show plots of the 
full SU(3) and projected potentials, for the case of one gauge-fixing/configuration. The 
potentials have been shifted by constants for convenient display, the significant quantity 
to compare being the slope of each potential at large R, which determines the string ten- 
sion. The string tension results obey the inequalities asu(3) > &u(i)xU(i) > <^z(3) > &mono- 
These inequalities hold at one gauge-fixing/configuration and become stronger when 
account is taken of gauge ambiguities. The most important result of this paper is 
that dsu(3) > °"c/(i)xc/(i)- When account is taken of gauge ambiguities at j3 = 6.0, the 
U{1) x U(l) string tension is smaller than the full SU(3) string tension by 20%, an 
amount well outside of error estimates. If a simulated annealing gauge-fixing algorithm 

|24| were used instead of the overrelaxation method used here, the discrepancy could 
be still larger. This situation for SU(3) is in contrast to SU(2), where the MAG and 
extrapolating to the global maximum of the functional (n g — > oo) produces a U(l) string 
tension very close to the full SU(2) result ||. 

It is conceivable that for SU(3), the discrepancy between <Jsu(3) an d cru(i)xU(i) could 
fade away with increasing /3 and become negligible in the continuum limit. This seems 
unlikely given that the ratio cfu{x)xU{i)I &su{3) is smaller at (3 = 6.0 than f3 = 5.9. Rather 
than moving into the region f3 > 6.0 using the Eq.(^|) form of the MAG, we intend to 
explore the possibility that the remedy lies in a generalized MAG such as Eq . fll5|) . The 
condition of Eq. (|15D seems to express in a natural way the idea that there is an effective 
adjoint Higgs field present, whose direction in the Cartan algebra may vary from point to 
point on the lattice. If it is accepted that the simple MAG form of Eq. (Q) cannot work for 
<S77(3), it is important to explore more general forms of the MAG. Failure to find any form 
of MAG which works would have serious implications for 't Hooft's proposal of abelian 
dominance of confinement for SU(N), N > 2. On the other hand, a quantitatively 
successful generalized MAG in SU (3) opens the possibility of a similar result for the 
indirect maximal center gauge, and a detailed study of the contribution to confinement 
of MAG monopoles and IMCG center vortices. These two degrees of freedom are clearly 
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intimately connected. 

With any generalized MAG, there will be a rather clear signal for improvement. An im- 
proved form of the MAG should produce a string tension larger than the full SU(3) result 
for one gauge-fixing/configuration. Gauge copies with higher values of the functional 
will inevitably lower the abelian £7(1) x U(l) string tension, hopefully to within error 
bars of the full SU(3) result. In contrast, a form of the MAG which for one gauge- 
fixing/configuration produces an abelian projected string tension smaller than the full 
SU(3) result is untenable. Gauge copies with higher values of the functional will only 
reduce it further. Our claim is that the MAG in the form of Eq.(§) is of this type. 

With regard to center projection, since Z{3) is a subgroup of U{1) x U(l), the group 
center should be accessible indirectly, i.e. after abelian projection. However, given the 
subtleties of the MAG for SU (3), a thorough study of the direct maximal center gauge for 
SU(3) is quite desirable. The DMCG avoids the question of what is the proper form of 
the MAG. We intend to study the direct maximal center gauge for SU(3) in future work. 
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Table 1: Full SU(3) results for string tension, Coulomb coefficient, and self energy at 
P = 5.90 and P = 6.0 



SU(3) 


P = 5.90 


P = 6.0 


Vol 


10 3 x 16 


16 4 


a 


0.068(3) 


0.050(1) 


a 


-0.36(2) 


-0.32(2) 


V 


0.68(2) 


0.66(2) 


R/a 


2 to 5 


2 to 7 



Table 2: MAG U(l) x U(l) results for string tension, Coulomb coefficient, and self energy 
at P — 5.90 and P = 6.0, one gauge copy/config. 



U(l) x U(l) 


P = 5.90 


P = 6.0 


Vol 


10 3 x 16 


16 4 


a 


0.063(3) 


0.045(2) 


a 


-0.09(2) 


-0.07(2) 


V 


0.20(2) 


0.19(2) 


R/a 


2 to 5 


2 to 7 
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Table 3: (3 — 6.0 MAG £7(1) x £7(1) string tensions for subsets of size n g out of 10 gauge 
copies 



n g 


Ct/(l)x[/(l) 


2 


0.042(1) 


4 


0.041(1) 


6 


0.041(1) 


8 


0.040(1) 


10 


0.040(1) 



Table 4: MAG monopole results for string tension, Coulomb coefficient, self energy, and 
current density, at f3 — 5.90 and f3 = 6.0, one gauge copy /configuration 



mono 


P = 5.90 


P = 6.0 


Vol 


10 3 x 16 


16 4 


a 


0.050(2) 


0.038(1) 


a 


0.03(2) 


0.045(1) 


V 


-0.04(2) 


-0.051(4) 


fm 


1.28(1) x 10" 2 


7.45(3) x 10" 3 


R/a 


2 to 6 


2 to 8 



Table 5: P — 6.0 MAG monopole string tensions and current densities for subsets of size 
n g out of 10 gauge copies 



n g 


&mono 


fm 


2 


0.036(1) 


7.25(4) x 10~ 3 


4 


0.035(1) 


7.08(4) x 10- 3 


6 


0.034(1) 


7.00(4) x 10~ 3 


8 


0.034(1) 


6.95(4) x 10^ 3 


10 


0.033(1) 


6.92(5) x 10- 3 
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Table 6: IMCG results for Z(3) string tension, Coulomb coefficient, self energy, and vortex 
density, at P — 5.90 and f3 = 6.0, one gauge copy/configuration 



IMCG 


P = 5.9 


P = 6.0 


Vol 


10 3 x f 6 


f6 4 


a 


0.060(3) 


0.040(2) 


a 


0.05(2) 


0.03(f) 


V 


-0.04(2) 


-0.02(f) 


V 


4.40(f) x 10~ 2 


3.01(1) x 10- 2 


R/a 


2 to 6 


2 to 8 



Table 7: P — 6.0 IMCG Z(3) string tensions and vortex densities for subsets of size n g 
out of fO fMCG gauge copies, best MAG gauge copy 



n g 


CZ(3) 


V 


2 


0.033(f) 


2.73(4) x f0~ 2 


4 


0.030(f) 


2.67(4) x f0~ 2 


6 


0.029(f) 


2.64(4) x f0~ 2 


8 


0.029(f) 


2.62(4) x fO" 2 


fO 


0.029(f) 


2.61(4) x 10~ 2 



Table 8: For cube dual to magnetic current link, r m {n) = fraction of time n faces are 
pierced by Z(3) flux 





p = 5.90 


P = 6.0 


r m (0) 


0.14(1) 


0.16(1) 


r m (2) 


0.74(1) 


0.74(1) 


r m (3) 


0.09(2) 


0.07(2) 


r m (4) 


0.02(2) 


0.02(2) 


r m (5) 


0.002(2) 


0.001(1) 


r m (6) 


0.001(1) 


0.001(1) 
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